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Abstract. We consider quadrature formulas based on interpolation using the ba- 
sis functions 1/(1 + t^x) (k = 1,2,3,...) on [—1,1], where are parameters on 
the interval ( — 1, 1). We investigate two types of quadratures: quadrature formulas 
of maximum accuracy which correctly integrate as many basis functions as possible 
(Gaussian quadrature), and quadrature formulas whose nodes are the zeros of the 
orthogonal functions obtained by orthogonalizing the system of basis functions (or- 
thogonal quadrature) . We show that both approaches involve orthogonal polynomials 
with modified (or varying) weights which depend on the number of quadrature nodes. 
The asymptotic distribution of the nodes is obtained as well as various interlacing 
properties and monotonicity results for the nodes. 



1. Introduction 
Suppose that we want to compute the integral 

■l 



f(x)w(x) dx 



l 



with w(x) a positive and integrable weight function on [—1, 1]; then one can use an 
interpolatory quadrature formula by approximating / using interpolation. In this 
paper we assume that w(x) > almost everywhere in [—1, 1] and use interpolation 
by means of the basis functions 

(i.i) 1 1 



1 + t\X 1 + t^x 1 + t n x 

where ti,t%, t$, . . . are parameters belonging to (—1, 1). We consider two approaches 
for choosing the interpolation nodes (quadrature nodes). The first consists of choos- 
ing the n nodes to be the zeros of the orthogonal function r n+ i obtained by orthog- 
onalizing the system (1.1) using the inner product 

(f,g) = / f{x)g(x)w(x) dx. 
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We refer to this as orthogonal quadrature. The second approach that we investigate 
is to choose the nodes in such a way that the quadrature formula is correct for as 
many basis functions as possible. This gives a quadrature with maximum accuracy, 
and we refer to it as Gaussian quadrature. It is well known that both approaches 
are the same when one is using interpolation by polynomials, but these choices 
result in different quadrature formulas when one uses interpolation with the basis 
functions (1.1). 

1.1 Orthogonal functions. It is well known that the zeros of the nth-degree 
orthogonal polynomial with respect to a positive weight function w(x) on [—1, 1] 
are all real, simple and in the interval (—1, 1), and that they separate the zeros 
of the orthogonal polynomial of degree n + 1 [26, p. 44-46]. These properties can 
be generalized to other orthogonal systems of functions {4>i : i = 1,2,...}, where 
orthogonality is with respect to the weight function w(x) on [—1, 1], provided they 
satisfy the Haar condition. The sequence of functions {fa : i = 1, . . . ,n} satisfies 
the Haar condition on [—1, 1] if and only if for every choice of n points xi G [—1,1] 
(i = 1,. . . ,n) with Xi ^ Xj whenever i ^ j one has deti<ij< n (fa(xj)) ^ (see, 
e.g., [14]). This condition is equivalent to saying that each linear combination 
4> = Yl7=i a i^i nas & t mos t n — 1 zeros in [—1, 1]. A system {4>i : i = 1, . . . , n} 
is an extended Haar system on [—1, 1] if and only if each linear combination <p = 
J27=i a i^i nas & t m °st n — 1 zeros in [—1, 1], counting multiplicities. Pinkus and 
Ziegler [25] and Videnskh [28] have investigated several properties of systems of 
orthogonal functions {4>i : i = 1, 2, . . . , n + 1} satisfying the Haar condition on 
[— 1, 1]. They have shown that 4> n +i has precisely n zeros in (—1, 1), the multiplicity 
of each zero is odd and the zeros of <p n and 4> n +i are strictly interlacing. 

Suppose that ti G (—1,1) (i = 1,2,...) and that ti ^ tj whenever i ^ j. 

Then the system j 1+ * ia; , 1+ t 2X , . . . , 1+t 1 +ix j is an extended Haar system on [—1, 1] 
because 

n+1 ^ ,\ 

oi-r) = ^ (i,o,[.r) 



nr=i 1 (i+^)' 

where p n (x) is a polynomial of degree at most n, and hence has at most n zeros 
on [—1,1], counting multiplicities. If we apply the Gram-Schmidt procedure to 

1 , i+ t2X 7 • • • | then we obtain a sequence of orthogonal functions on [—1, 1] 

with respect to the weight function w(x). We denote these orthogonal functions by 

1 n c- 

r„+i(*i, . . . ,t n+1 ;x) = — — h Y] \ , 

l + t n+1 x ^1 + Ux 

and, as mentioned before, this rational function has precisely n zeros in (—1, 1). It 
is well known that 

1 n a- 2 

> — w(x) dx 

l + t n+1 x f-^l + Ux 



mm 

(oi,... ,a„)6l 



is obtained when 

1 n 
1+^-^ = ^ 

so that r n+ i(ti, ... , t n ;x) is the error function for the best L2(w)-approximation 

nf hv linpar rnmhina.tinns nf / -. — - — - — - — - — - — X 
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1.2 Gaussian quadrature. Suppose the system of functions {4>i : i = 1,2, . . .} 
is an extended Markov system on [—1, 1], i.e. : i = 1, 2, . . . , k} is an extended 
Haar system on [—1, 1] for each k G N. If one considers a quadrature formula with n 
nodes which is of maximum accuracy, then one wants to integrate the functions (pi 
(1 < % < m) correctly for m as large as possible. If 4>i{x) = x % ~ x , then m turns out 
to be 2n, and the nodes are the zeros of the orthogonal polynomial of degree n for 
the weight function w(x) (Gauss-Jacobi quadrature). Krein [15] showed that when 
{<pi : i = 1, 2, . . . , 2n} is a Haar system on [—1, 1], then there exists a unique set of 
n nodes Xi G (—1,1) and n strictly positive weights fa such that for k = 1, 2, . . . , 2n 



and this gives maximum accuracy. Karlin and Pinkus [12] have generalized this by 
allowing nodes with multiplicities and there are further generalizations by Barrow 
[2] and Bojanov, Braess and Dyn [3]. 

The system of rational functions (1.1) is an extended Markov system on [—1, 1], 
hence the properties mentioned above are all valid. In particular it follows that the 
nodes for Gaussian quadrature are all in (—1, 1) and the corresponding quadrature 
weights are all strictly positive. 



In this section we will make the simple observation that the nodes and weights 
for both, orthogonal quadrature and Gaussian quadrature, based on the rational 
functions (1.1) are closely related to the zeros and Christoffel numbers for orthog- 
onal polynomials on [—1, 1] with a weight function that depends on the number of 
nodes. This observation is not new (see e.g., Lopez and Illan [18] [19] for Gaussian 
quadrature), but this interaction makes it possible to use results from the theory 
of orthogonal polynomials to obtain useful properties of the nodes and weights for 
quadrature based on rational interpolation. 

Theorem 1. Let {4>k} be the set of basis functions (1.1). 

(1) Suppose all the ti are different; then the n nodes for orthogonal quadrature 
are the zeros Xj jn (1 < j < n) of the orthogonal polynomial p n (x) with 
respect to the weight function 



(1.2) 




i=l 



2. Connection with orthogonal polynomials 



(2.1) 



w n (x) 



w(x) 



-1 < X < 1, 



n n (x)n n+1 (x) ' 



where 



n 



(2.2) 



n 



(x) = Y[(i+ux). 



If t\ =0 then (1.2) holds for k = 1, 2, . . . , n + 1 with Xj = Xj tU and 



Pj(ti, . . . , t n+ i) 



where A J;n are the Christoffel numbers corresponding to the orthogonal poly- 
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(2) Suppose all the t{ are different; then the nodes for Gaussian quadrature are 
the zeros Xj jU (1 < j < n) of the orthogonal polynomial p n (x) with respect 
to the weight function 

(2.3) w n (x) = , -1< x < 1, 

and (1.2) holds for k = 1, 2, . . . , 2n with Xj = Xj jn and 

Pj(tl, ■ ■ ■ ,^2n) = K2n{Xj,n)^j,m 

where A J;n are the Christoffel numbers corresponding to the orthogonal poly- 
nomial p n (x) for the weight function (2.3). 

Proof. 

(1) The orthogonal functions obtained by orthogonalizing the system of basis 
functions (1.1) on [—1,1] with weight function w(x) are of the form 



n+l 



r n+1 (tu • • • , Wi; x) = ^ Y+t^ = ^{x) {an+1 = 1} ' 

where p n (x) is a polynomial of degree n. The orthogonality implies that 

11 Pn{x) 



f P , V — ^ — w(x) dx = 

J_ x 7T n+1 (x) 1 + tjX 



for every choice of parameters ci, c 2 , . . . , c n . This is the same as saying that 

1 w(x) 
p n (x)q n -i(x) — -—dx = 



n n (x)n n+1 (x) 

for every polynomial q n -i(x) of degree at most n — 1. This shows that p n (x) is the 
orthogonal polynomial of degree n for the weight function (2.1). The Gauss-Jacobi 
quadrature for the weight function (2.1) gives 

a a\ f 1 q 2n -i(x)w(x) <A 

(2 ' 4) / M*K + iM dx = g A -*»-' (l -»' 

where the zeros of p n (x) and q2n-i( x ) is a polynomial of degree at most 

2n — 1. If ti = then 7r n (x) is a polynomial of degree n — 1. In this case choose 
<?2n-i(£) = 7r n (x)7r n+ i(a;)/(l + t fc x) to obtain 



/_ 



W ( X ) dx = Y" X J,n^n(x Jjn )TV n+1 (x Jjn ) k=12 n + l 

1 -I- 1 1. T -• _ 



1 1 + ^fc^ . =1 1 + tkXj jn 
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(2) Let Xi = Xi(ti, . . . , t 2 n) and A = (3i{t u . . . , t 2 n) (i = 1, 2, . . . , n). In order 
to obtain the nodes Xi and weights one needs to solve the system of equations 

1 T^T-^ = Et^-' j = l,2,...,2n. 
_i 1 + f-f 1 + t 7 -Xi 

This is equivalent to finding Xj and (i = 1, 2, . . . , n) such that 

J-l K2n{x) 7T2n{Xi) 

for every polynomial qi n -\ of degree at most 2n — 1. From the theory of orthogonal 
polynomials (Gauss-Jacobi quadrature) we then know that the nodes Xi are the 
zeros of the orthogonal polynomial of degree n for the weight function (2.3) and 
then 

2n ~* w(x) cte = Ai ) „g 2 „_i(x < ). 

Choose g 2 n-i(^) = 7T2n(^)/(l + *fc^) (fc = 1, . . . , 2n) to obtain the assertion. □ 
In case the weight function w(x) is equal to either 

1 

w(x) = \ Vl - x 2 , 

1-x 
1+£C' 

the orthogonal polynomials for the modified weight functions w n (x) given by (2.1) 
or (2.3) can be found explicitly and are known as Bernstein-Szegd polynomials [26, 
P- 31]. 

The nodes of both, orthogonal quadrature and Gaussian quadrature, are equal 
to the zeros of the orthogonal polynomial p n (x) with respect to a weight function 
depending on n. Therefore, the distribution of the quadrature nodes for orthogonal 
quadrature and Gaussian quadrature is given by the distribution of the zeros of 
orthogonal polynomials with a weight function depending on the degree n. The 
asymptotic zero distribution and various asymptotic results for such orthogonal 
polynomials have been investigated recently by Mhaskar and Saff [22]-[24], Gon- 
char, Lopez and Rakhmanov [11] [21] and Lopez [16] [17] [20]. An almost immediate 
corollary to their results is: 

Theorem 2. Suppose that the asymptotic distribution of the parameters {U : i < 
n} is given by a measure v on [—1,1], i.e., for every continuous function f on 
[—1,1] we have 

(2.6) lim -J2f(tj)= / f(x)dv(x). 

n^oo n ' J_i 
J = l 

If ' v = pd-i + qdi + rdo + [1 — (p + q + t)]vq, where 5 a is the unit measure with all 
its mass concentrated at a, and p,q,r>0,p + q + r<l and 



(2.7) / log It] du (t) < oo, 



r 
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then the asymptotic distribution of the nodes for both, orthogonal and Gaussian 
quadrature, is given by the measure 

(2.8) fj, = p6 1 + q5-i + rfj, + [l - (p + q + r)]fi b , 

where /j,q is the arcsin measure on [— 1, 1] with weight function 

H {x) 



7T VT - X 2 

and (ib is an absolutely continuous measure on [— 1, 1] with weight function 

(2.9) /li(l) = I_i_£^^ (t) . 

Proof. For orthogonal quadrature the weight function for the orthogonal polyno- 
mials is w n (x) = w(x)e~ 2niPn( - x ^ with 

<Pn(x) = ^ [log7T n+ i(x) + log7T n (x)] , 

whereas for Gaussian quadrature it is w n (x) = w{x)e~ 2n ^ n ^ with 

V'n = ^-^ogn 2n (x). 
2n 

By the weak convergence (2.6) we have 

(2.10) lim (p n (x) = lim ip n (x) 

= p log(l — x) + q log(l + x) + [1 — (p + q + r)] J log(l + tx) dv (t), 

uniformly on closed sets of (—1,1). By [11, Theorem on p. 124], the asymptotic 
distribution of the zeros of the corresponding orthogonal polynomials is given by a 
measure fx on [—1,1] which is the unique solution of the integral equation 

f 1 1 

(2.11) / log, -dfi(t) + <p(x) = C, xe[-l,l], 

J-i \x — t\ 

where ip(x) is the right-hand side of (2.10) and C is a constant. We easily find 

(p(x) =plog(l-x) +qlog(l + x) - [1 - {p + q + r)] J log j_ dv (-l/t) 

+ [1 - (p + q + r)} J \og\t\dv (t). 

The last term in this expression is a finite constant by (2.7) and may be absorbed 
by the constant in (2.11). Notice that (p(x) may become — oo near ±1, which at 
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of the proof in [11] shows that the case under consideration is also covered. The 
integral 

(p b (x)= log | - dv (-l/t) 

</(-oo,-l)U(l,oo) \X — t\ 

is the logarithmic potential of the mass distribution vq{— 1/t) outside [—1,1]. The 
integral equation 



1 1 



log _ d/i b (t) - ip b (x) = C, £€[-1,1], 

has a unique solution fib which is the 'balayage' of the mass distribution vq{— 1/t 2 ) 
outside [—1, 1] to the interval [—1, 1] (see, e.g., [21, pp. 126-127]), and this measure 
fib is absolutely continuous with weight function given by (2.9). The arcsin distri- 
bution has a constant logarithmic potential on [—1, 1], which can be absorbed into 
the constant C in (2.11). Hence, the unique solution of (2.11) is the measure given 
by (2.8). □ 

Let us take a closer look at the limiting zero distribution measure \i given by 
(2.8). If the measure v has all its mass at the point a G (—1, 1) then 



(2.12) fj?(x) 



l VT 



a 



2 



7T VT - x 2 l + ax 



If a = 0, then the limiting zero distribution is the well-known arcsin distribution. 
For a = — 1 the limiting zero distribution is a degenerate distribution at 1 wheras 
for a = 1 the limiting zero distribution is concentrated at —1. If a is close to 
1 in (2.12), then // has even more mass concentrated around —1 than the arcsin 
distribution. The measures with weight functions given by (2.9) or (2.12) are thus 
somewhere between an arcsin measure and degenerate measures at ±1. 

Even though the asymptotic distribution of the nodes for orthogonal quadrature 
and Gaussian quadrature is the same, we see that the n nodes for the Gaussian 
quadrature formula depend on 2n parameters U whereas the n nodes for orthogonal 
quadrature depend on n + 1 parameters ti. The nodes for Gaussian quadrature 
therefore use more (asymptotic) information of the basis functions. 

3. Properties of the nodes 

The orthogonal functions r n (ti, ... ,t n ;x) (n = 1,2,...) are an extended Haar 
system on [—1, 1], and therefore the properties of the zeros given by Pinkus, Ziegler 
and Videnskh are valid. For this special Haar system additional useful properties 
can be proved: 

Theorem 3a. 

(1) The zeros of r n (tx, . . . ,t n ; x) and r n+ i(£i, . . . , t n , t n+1 ; x) on [-1,1] are 
strictly interlacing. 

(2) The zeros of r n {t u ... ,t n _ u t n - : x) and r n (ti, . . . , £ n _i, t n+1 ; x) on [-1, 1] 
are strictly interlacing. 

(3) The zeros of r n {t u . . . , t i+1 , . . . , t n , tf, x) (1 < i < n) and r n {t u . . . , 
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Proof. Property (1) is proved in Pinkus and Ziegler [25] and Videnskri [28], and 
property (2) can be found in Pinkus and Ziegler [25]. For property (3) we notice 
that 

T-n (tl , • • • , t n , X) T n — \ (tl , . . • , t i—\ , 5 • • • jt n , X) 

KlTn— 1 (tl , • • • , ti— l , , . . . , t n —\ , tjj 2?) 

and 

(tl , • • • i ti — 1 ? ) • • • 5 tn, tj, 2-) ?"n — 1 (tl ? • • • > ^i— 1 > > • • • , tn— 1 , ti, x) 

~ ^2^n-l (tl, • • • , ti-1, ti+1, • • • j t n -l, t n 'i x )i 



where 



^'"n— 1 (tl , • • • , ti— 1 , ti+1 1 • • • i t n ? *^) , 1+i, 



«i = 7" : : : — , , |2 



and 



r n— l(tli • • • , tj_ 1? ti+li ■ ■ ■ ,t n _ i,ti,x) 

[ fn— 1 (tl j • • • , ti— l, ti+l, • • • , t n — 1 , ^i! ^0 , T+t~x ) 

^ 2 = ii '( : ^ T ~ 1 

l^n— 1 (tl 1 • • • 5 ti— 1) ti+lj • • • jt n ,X 



12 

Let Xi (i = 1, 2, . . . , n — 1) be the zeros of r n (ti, ... , t n ; x); then 

(3.1) r n (t\, . . . , ti— i, ti+i, • • • it n ,ti,Xj^) 

= r n -i(t\, ... , ti— i, ti+i, • • • , t n _i, ^; — «i«2)- 

If we compute the L 2 -norm of r n (ti, . . . ,t n ;x), then 

||^n(ti) • • • , t n ; x) H2 = I r n (ti, . . . , t n ; x), - — I 

V l + t n xj 

/ 1 

I Tn— 1 (tl , . . . , ti — 1 , ti_|_i , ... , t n , X) , 



1 + t n £, 

1 



' 1 + t n X 



K>1 y^n— l(tl>... 5 ti— 1) ti+lj • • • ) t n _ i,ti;x) 

= ||r n _i(ti, . . . ,tj_i,t i+ i, . . . ,t n ; x) - Avik 2 ) 

so that 1 — > 0. Then property (3) follows immediately from (3.1) and property 
(1). □ 

We have similar results for the nodes corresponding to Gaussian quadrature. 
Now the order of the parameters ti, . . . , t 2n is irrelevant. 

Theorem 3b. Let Xi(ti, . . . , t 2n ) (i = 1,2, ... ,71) 6e t/ie quadrature nodes for 
maximum accuracy. 

(1) The n nodes Xi(ti, . . . , t 2n ) and the n + 1 nodes Xj(ti, . . . , t 2n +i, t 2n + 2 ) are 
strictly interlacing. 

(2) TTie n nodes Xi(ti, . . . , t 2n -i, t 2n ) and t/ie n nodes Xi(ti, . . . , t 2n -i, t 2n +i) 
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Proof of (1). Let x Q {t u . . . ,t 2n ) = -1 and z n+ i(ti, . . . ,t 2n ) = 1- 
Suppose that the nodes Xi(t\, . . . , t 2n ) and xi(ti, . . . , ^2n+2) are not strictly inter- 
lacing; then either 

(a) Xi(ti, . . . ,t2n) = Xj(ti,... ,^2n+2) for some i and j with 1 < i < n and 
1 < j < n + 1 or 

(b) Xj(ti, . . . ,t 2n +2), Xj+l(h, . . . ,t 2 n+2) £ ((Xi(tl, . . . ,t2n),Xi+l(tl, ■ ■ ■ ,*2n)) for 

some i and j with < i < n and 1 < j < n. 

We will show that both (a) and (b) are impossible. 

Case (a): Suppose that there are i nodes Xi(ti, . . . , t2n+2) 5 i £ {ii, ^2, • • • , u}, such 
that Xi(ti, . . . ,t2n+2) £ {xj(ti, . . . , t2n) '■ 1 < 3 ' < then obviously £ < n, and 
thus there exists a /c such that Xfc(ti, . . . , ^2n+2) ^ {^'(^1? • • • ? ^2n) : 1 < J ' < 
Let be a linear combination of the first 2n basis functions given in (1.1) such 
that (xi(ti, . . . , tin)) = for i = 1, 2, . . . , n and . . . , ^2n+2)) = for 

z = 1,2,... , 2n+ 1 but i ^ k. Then we already have 2n — £ equations to determine 
0. In addition we add the equations <f>(zi) = for i = 1, . . . , £ — 1, where the ^ 
differ from the 2n — £ points that we already used. For such we have 

/l n 
4>{x)w(x) dx = ^ @i(ti, . . . ,t2n)<l>(Xi(tl, ■ ■ ■ ,t 2n )) = 0. 
_1 i=l 

But we also have 

/l n+l 
(j)(x)w(x) dx = J~]/3j(ti, . . . ,t 2n +2)(p(Xi(ti, . . . ,*2n+2)) 
- 1 i=l 

= fik(tl, . . . , t2n+2)4> (Xk(tl, • • • , ^2n+2)) 7^ °, 

because can have at most 2n — 1 zeros and (3k > 0. This gives a contradiction. 
Case (b): Let be the linear combination of the first 2n basis functions given in 
(1.1) such that (xi(ti, . . . , t 2n )) = for % = 1, 2, . . . , n, and . . . , ^2n+2)) = 

for z = 1, 2, . . . , n + 1 but with i ^ j and i ^ j + 1. This gives 2n — 1 conditions 
to determine 0, and for such <f> we have 

0(x)«;(a;) da; = • • • > ^2n)</> (^(^i, • • • , *2n)) = 0. 

On the other hand we also have 

/l n+l 
<f>(x)w(x) dx = ^2 A(*i> • • • > ^2n+2)0 (a;i(*i, . . . , t 2n + 2 )) 
_1 i=i 

= /3j(tl, . . . ,t 2 n+2)<P(Xj(ti, . . . ,t 2n +2)) 

+ (3j + i(ti, . . . , t2n+2)<j>(%j+l(tl, ■ ■ ■ ,t2n+2)) ■ 

By assumption we have sign (xj(ti, . . . ,t 2n+2 )) = sign . . . ,t 2 n+2)), 

and since the weights ^ and are positive, we have a contradiction. 

Proo/ of (2). Suppose that Xi = x^tx, . . . ,t 2n ) and y t = Xi{t u . . . Mn-iMn+i) 
are not strictly interlacing; then either 
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(b) Xi = yj for some i and j with 1 < i < n and 1 < j < n; or 

(c) yj, yj+i G (x^ Xi+i) for some i and j with < z < n and 1 < j < n — 1. 
Case (a): The system of basis functions (1.1) is a Haar system on [—1, 1], hence we 
can find a function fk in the linear space spannen by the first 2n — 1 basis function 
in (1.1) such that 

fk(xi) = $i,k, i = l,2, ... ,n, 
/*(&) = 0, z = l,2,...,n-l 

where £j G (—1, 1) \ {x\, x 2 , • • • , x n } are n — 1 arbitrary points. Then on one hand 

,1 n 

f k (x)w(x) dx = ^2 A(*i> • • • » hn)fk{Xi) = Pk(h, . . . , £ 2 n), 
and on the other hand (since = yi for 1 < i < n) 

/l n 
f k (x)w(x)dx = y~]/3j(ti, . . . ,t 2 n-l,t 2 n+l)fk(Xi) = Pk(tl, . . . , *2n-l , fon+1 ) ■ 

Hence (3i{t\,... Mn-iMn+\) = ■ ■ ■ Mn) for z = 1, 2, . . . , n. Therefore the 

quadrature formula will give the exact result for every function with % = 

1,2,... , 2n + 1. Let be a linear combination of these 2n + 1 basis functions such 
that = and 4>'(xi) = (i = 1,2, . . . ,n); then does not changes sign on 

[—1,1] because <f>(x) = ±p^(x) /7V2 n +i(x) with p n a polynomial of degree n and 
7r 2n +i given by (2.2). This implies that 

(j)(x)w(x) dx 7^ 0. 

i-i 

This contradicts 

/! n 
4>{x)w{x) dx = (3i<j)(xi) = 0. 
_1 i=i 

Cases (b) and (c): one can obtain a contradiction in a way similar to cases (a) 
and (b) of (1). □ 

The following result of A. A. Markov (see, e.g., [26, Thm. 6.12.1 on p. 115]) allows 
us to give more information on the behavior of the nodes Xi(t\, ... ,t n ,t) for or- 
thogonal quadrature, and xi{t\, . . . , ti n -\, t) for Gaussian quadrature, as functions 
of t G (-1,1): 

Theorem. (Markov) Let w(x, t) be a weight function on [a, b] depending on a pa- 
rameter t such that w(x, t) is positive and continuous for a < x < b and c < t < d. 
Suppose that w t (x,t) = ^w(x,t) exists and is continuous for a < x < b and 
c < t < d, and that the moments 

rb 

/ x h w t (x,t) dx, k = 0, 1, 2, . . . ,2n — 1, 

J a 

converge uniformly in every closed subinterval c' < t < d' of (c,d). If the zeros of 
the orthogonal polynomial p n (x,t) with weight function w(x,t) are denoted by Xi(t) 
(i = 1,2, ... ,n) and if w t (x,t)/w(x,t) is a strictly increasing (decreasing) func- 
tion of x & (a,b), then Xi(t) is a continuously differentiate and strictly increasing 
(decreasing) function of t G (c, d) for every fixed i. 
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Theorem 4. 

(1) If Xi(ti, . . . ,t n ,t) is the ith node for orthogonal quadrature with n nodes, 
then Xi(ti, . . . ,t n ,t) is a continuously differentiate and strictly decreasing 
function of t G (—1,1). 

(2) Every node Xi(ti, . . . ,t2 n -i,t) for Gaussian quadrature with n nodes is a 
continuously differentiate and strictly decreasing function of t G (—1,1). 

Proof. By Theorem 1 (1) we know that Xi(ti, ... ,t n ,t) is the ith zero of the or- 
thogonal polynomial p n (x, i) of degree n for the weight function 

/ \ w(x) 
w(x,t) = 



nl(x)(l + tx)- 
A simple computation gives 

Wt(x,t) x 
w(x, t) 1 + tx 

and this is a decreasing function of a; G (—1,1). The result thus follows by Markov's 
theorem. 

Similarly, by Theorem 1 (2), for the nodes Xi(ti, . . . , tm-u t) of Gaussian quad- 
rature one can use Markov's theorem for the weight function 

/ \ w(x) 
w(x,t) = 



7r 2 „_i(x)(l + tx)' 

to obtain the desired result. □ 

We can also give some monotonicity results for the first and last quadrature 
coefficients. Let /%(ii,... ,t n ,t) be the quadrature coefficients corresponding to 
the nodes Xi(ti, ... ,t n ,t) for orthogonal quadrature with n nodes, and let t\ = 0; 
then the first quadrature coefficient fii(ti, . . . ,t n ,t) is a strictly decreasing con- 
tinuously differentiable function of t e (—1,1), and the last quadrature coeffi- 
cient P n (ti, ■ ■ ■ ,t n ,t) is a strictly increasing continuously differentiable function 
of t G (—1,1). Similarly, if . . . , tm-\i t) are the quadrature coefficients 

for Gaussian quadrature with n nodes (with t\ = 0); then fii(ti, . . . ,t2 n -i,t) 
is a strictly decreasing continuously differentiable function of t G (—1,1) and 
(3 n (ti, . . . ,t2 n -i,t) is a strictly increasing continuously differentiable function of 
t G ( — 1, 1). The proofs of these results are rather technical (see [27]) and since only 
the extreme quadrature weights are considered we decided not to include them here. 

4. Convergence results 

If the integral 

l 

f(x)w(x) dx 

is approximated by the orthogonal quadrature sum with n nodes, 



n 

pijti, ■ ■ ■ ,t n+1 )f{xi{t u ... ,t n+1 )), 
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we denote the error by 



/l n 
f(x)w(x)dx-^2p i (t 1 ,... ,t n+1 )f(x i (t 1 ,... ,t n+1 )). 

_1 3 = 1 

The error of the Gaussian quadrature formula is denoted by 

/l n 
f(x)w(x) dx-^Pifa,... ,t 2n )f(x i (t 1 ,... ,t 2n )). 

- 1 3=1 

It is clear that these quadrature formulas are only going to be relevant if for a 
large enough class of functions / one has E°f — > as n — > oo and/or E®f — > 
as n — > oo. It is well known that the sequence of basis functions (1.1) is dense in 
C[— 1, 1] with respect to the uniform norm if and only if 

oo 

(4.1) - |cfc|) = oo, 



k=l 



where Ck = — 1/tfc — a/V^1 — 1 [1; P- 254, §7] 1 . The following result is straightfor- 
ward, but we mention it for the sake of completeness. 

Theorem 5. Let t±, t 2 , £3, . . . be distinct parameters in (—1, 1) and <f> n +i(x) be the 
best approximation of f E C[— 1, 1] using linear combinations of , i+t 2 x ' • • • > 



(1) = 0, then for orthogonal quadrature, there exists a positive constant M 
not depending on n such that 

K/|<M||/-0 n+1 || oo . 

As a consequence, E°f — > as n — > oo /or every bounded Riemann inte- 
grate function on [—1, 1] whenever (4.1) holds. 

(2) For Gaussian quadrature we have similarly 

\E°J\ < M\\f - fonWoo- 

As a consequence, E^f — > as n — > oo /or every bounded Riemann inte- 
grate function on [—1, 1] whenever (4.1) holds. 

Proof. The result for orthogonal quadrature can be proved as in the case of poly- 
nomial interpolation [8, pp. 126-129]. For Gaussian quadrature we only need to 
show that Yl7=i A(*i> • • • ? ^2n) is bounded. If £i = 0, then 



n „i 

• • »*2n) = / V)(x)dx. 



1 There is a misprint in this reference: 'interior' should be replaced by 'exterior'. We thank P. 
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If £1 ^ 0, then 



/" w(x) = A 

7_1 1 + ^ l + ti£ 

The function 1/(1 + tix) has the minimum value 1/(1 + |ti|) on [—1, 1], hence 

1 w(x) 



1 f 



+ t\X 



which shows the desired boundedness. □ 



Let / G C[— 1, 1] and denote by L n f the rational interpolant of / at the points 
Xi(t!,... , t n+1 ) (i = 1, 2, . . . , n) using the basis functions j^^, 1+7^, • • • , 
The previous theorem implies that 



lim / (L n f)(x)w(x) dx = / f(x)w(x)dx, 

whenever (4.1) holds. The rate of convergence depends on how well / can be 
approximated in the uniform norm by the basis functions (1.1), and this depends on 
the choice of the parameters ti. We now show that rational interpolation converges 
in L 2 (w) and thus also in Li(w), thereby extending the famous Erdos-Turan result 
for polynomial interpolation [6]. 

Theorem 6. Suppose t\ = and t 2 , ^3, • • • are distinct parameters in (—1, 1) such 
that (4.1) holds. If f G C[— 1,1] and if L n f(x) = Yl7=i i+t-x interpolates f at the 
zeros Xi(ti, . . . ,t n+1 ) o/r n+ i(*i, . . . ,t n+1 ;x), then \\L n f - f\\ 2 -> as n -> oo. 

Proof. Let <f> n (x) be the best uniform approximation of / using the basis functions 

{i+b'--- 'T+fcJ; then 

\\L n f - fh < 11/ - 4>nh + \\L n f - 4>nh- 

Clearly, ||/ — 4> n \\2 — *• as n — > oo. Now [(L n f)(x) — <p n (x)] 2 is a linear combination 
of the functions x^ /[tt 2 (x)] (j = 0, 1, . . . ,2n — 2), and by (2.4) the quadrature 
formula is therefore exact for [(L n f)(x) — <p n (x)] 2 . This gives 



\L n f ~ (pn\\l = ^2Pi[(L n f)(xi) - <fi n (Xi)f 
i=i 

n 

= X>[/(zi)-*»(*i)] 2 

i=l 

< \\f -KWlo J ^w{x)dx. 



Thf rpsnlt nnw fnllnws hpr.a.nsp II f — r/>„ II — > as r?. — ► no. i 
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5. Generalizations 



All the properties and results for orthogonal quadrature that we have given above 
are also true for the orthogonal system obtained by orthogonalizing the system of 
basis functions 



(5.1) 







1 + Ux ' dt 1 + tx 



-1 



t=u 



<9£Mi-l l + tx 



i — 1, . . . , k, 



t=u 



where k G N, \ii G N and X)i=i A 4 * = n + 1- If the parameters ^ are all distinct and 
in (—1,1), then this system of functions is an extended Haar system on [—1,1]. If 
t {£i, ... ,tk} then 



k m — I 



i=l j=Q 



t=U 



is orthogonal to the system (5.1) with respect to the weight function w(x); ift = tg 
for some £ with 1 < I < k then 



_ (Am) AfJ-e-i) Avt+i) AfJ-k) \ 



dt^ 1 + tx 



k — 1 
t—te i—i j—o 



1 



1 dP 1+tx 



t=u 



is orthogonal to the system (5.1) with respect to the weight function w(x). Because 
of the Haar property it follows that these orthogonal functions r n+ i have precisely 
n zeros in [—1, 1] which are all simple [25], [28]. These orthogonal functions share 
the same properties as the orthogonal functions for the system (1.1): the reason for 
this is that Theorem 1 can be extended to this system. Similarly all the results and 
properties for Gaussian quadrature that we have given are also true for the system 
of basis functions (5.1) with k G N, \Xi G N and X)i=i A*i = ^n. The only things that 
change are the weight functions (2.1) and (2.3), which need to be appropriately 
modified. 



Theorem 7. 

(1) The zeros of r n+1 (^ Ml) , . . . , t ( ^ k) ; 

polynomial p n (x) of degree n for the weight function 



tT K/ ;x) are also the zeros of the orthogonal 



w(x) 



(2) The n nodes for Gaussian quadrature with the system (5.1) are the zeros of 
the orthogonal polynomial p n (x) of degree n for the weight function 

w(x) 



n?=i(i+*i*) w 



This theorem allows us to use Markov's theorem to deduce properties of the 
nodes and the weights for orthogonal quadrature, but since everything is along the 
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6. Numerical examples 

Prom Theorem 5 it follows that rational quadrature rules will be effective when- 
ever the function to be integrated can be approximated well using the basis func- 
tions (1.1). Gautschi [9] has shown that functions which have an infinite number of 
poles but are regular otherwise, can be integrated most effectively by Gauss rational 
quadrature rules by choosing the parameters ti in an optimal way. Gautschi [10] ex- 
plicitly used such quadrature rules for the computation of generalized Fermi-Dirac 
and Bose-Einstein integrals. 

For functions with most mass concentrated near the endpoints of the interval 
[—1,1], one could choose the parameters ti in such a way that more quadrature 
nodes are situated near ±1. If we take the parameters 

(6.1) t i = i_J_ ten, 

then there will be more quadrature nodes near —1 than near +1. In fact, from The- 
orem 2 it follows that the asymptotic distribution of the nodes for both, orthogonal 
quadrature and Gaussian quadrature, is degenerate at the point —1. Let GR(n) be 
the n-point Gaussian quadrature rule, and OR(n) the n-point orthogonal quadra- 
ture rule. The nodes and weights for n = 6 are given in Table 1. Our computations 
have been made on an IBM 3090/400E VF using Brent's multiple precision package 
[5], and an independent computation was made using mathematica™ [29] 2 . 



nodes for GR(6) 


weights for GR(6) 


-0.9797390942708352 
-0.8853794251591486 
-0.6822351336410264 
-0.3156675377072605 
0.2408527285476740 
0.8155273184304977 


0.0528758827013522 
0.1411615118844550 
0.2748067575758422 
0.4657849717765712 
0.6221630733806293 
0.4432078026811501 


nodes for OR(6) 


weights for OR(6) 


-0.9736320979338328 
-0.8537169072027923 
-0.6094091127142633 
-0.2057016948376719 
0.3414560761423378 
0.8474273771128526 


0.0685126325838336 
0.1760476819554412 
0.3192517203251832 
0.4878639628808742 
0.5765658940369015 
0.3717581082177663 



Table 1: nodes and weights for the quadrature rules GR(6) and OR(6) 

The nodes and weights for GR(n) and OR(n) have been computed by using the 
zeros and Christoffel numbers for the orthogonal polynomial of degree n with the 
corresponding weight (2.1) for orthogonal quadrature and (2.3) for Gaussian quad- 
rature. These orthogonal polynomials have been computed using their three-term 
recurrence relation. The recurrence coefficients were obtained using Chebyshev's 
algorithm with modified moments and the latter were computed recursively by 
decomposing the rational weight function into simple fractions. 
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Consider the integral 




dx, u>0. 



J-i 



This function has more mass near —1 than near +1, especially for large u. In 
Table 2 we have given the absolute relative errors for the Gaussian quadrature 
sum and the orthogonal quadrature sum corresponding to the rational functions 
with parameters given by (6.1), together with the Gauss-Legendre quadrature sum 
GL(n) with n nodes. The rational quadrature rules behave in a similar way and 
give better results than the Gauss-Legendre quadrature rule. 



u n 


GR(n) OR(n) GL(n) 


5 2 
6 

10 
14 


0.664(-l) 0.121 0.394 
0.261(-5) 0.207(-6) 0.161(-4) 
0.215(-8) 0.842(-10) 0.212(-5) 
0.494(-ll) 0.514(-13) 0.112(-5) 


25 2 
6 

10 
14 


0.978 0.988 0.999 
0.118(-4) 0.255(-4) 0.206 
0.397(-8) 0.467(-8) 0.201(-2) 
0.206(-12) 0.122(-12) 0.153(-4) 



Table 2: relative errors for Ii 



The rational quadrature rules perform well because the system of rational func- 
tions with parameters (6.1) emphasizes the behaviour of the integrand near —1, and 
even though the quadrature nodes have a degenerate asymptotic distribution at —1, 
the parameters still satisfy (4.1) so that the basis functions are dense in C[— 1, 1]. 
The fact that the poles of the rational functions are all close to —1 also explains why 
both rational quadrature rules give similar results: adding extra rational functions 
is not giving a lot of extra information about the integrand. 

As a second example we consider the integral 

h = I' dX 

J-i y/(x + 3)(x + 2)' 

The integrand has a branch cut on the interval [—3, —2] and is in fact a Stieltjes 
function given by 

i _ i r 2 _L dt 

v /(x + 3)(x + 2) ~ 7v J_ 3 x-t v /( 3 + t)(_2-t)' 

Stieltjes functions can be well approximated by rational functions with poles on 
the cut. Therefore we will choose the parameters ti in such a way that the poles of 
the basis functions (1.1) are distributed nicely on [—3, —2]. In order to do this we 
observe that the zeros of a Chebyshev polynomial of the first kind T 3 k (x) are also 
zeros of Chebyshev polynomials T% m (x) for every m > k. Hence we take t\ = (so 
that Theorem 1 can be applied) and t^-, £3, • • • , £3™+! equal to 



-2 
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where 0^3™ are the zeros of the Chebyshev polynomial T 3 ™ (x) ordered in such a way 
that xi^m, £2,3™, #3,3™ are the zeros of Ts(x) in decreasing order; #4,3™, . . . , 0:9,3™ 
are those zeros of Tg(x) which are not zeros of Ts(x), again in decreasing order; and 
recursively we let be the zeros of T 3 k (x) which are not zeros 

of T 3 fe-i(x), in decreasing order. In this way the poles of the basis functions (1.1) 
are dense on [—3, —2] and they are distributed on [—3, —2] according to the arcsin 
measure on this interval. In Table 3 we have given the relative errors for the two 
rational quadrature rules and for the Gauss-Legendre rule. 



n 


GR(n) 


OR(ra) 


GL(ra) 


2 


0.587(-6) 


0.168(-4) 


0.366(-2) 


6 


0.161(-19) 


0.262(-16) 


0.567(-7) 


10 


0.986(-31) 


0.132(-29) 


0.116(-11) 



Table 3: relative errors for I2 



It is clear that the rational quadrature rules give a much better approximation to 
the integral I2 than the Gauss-Legendre rule. The reason is obvious: the integrand 
is a Stieltjes function which can be approximated much better by rational functions 
than by polynomials, especially when we choose the poles of the rational functions 
on the interval [—3, —2] on which the Stieltjes function has its singularities. Also 
the Gaussian quadrature rule is slightly better than the orthogonal rule, because 
the Gaussian rule with n nodes already uses 2n — 1 poles on the interval [—3, —2], 
whereas the orthogonal rule only uses n poles. 

Finally we consider the integrals 

h= . '. r eft, h= — 

J_ 1 sin(7!r/a;J J_ 1 \sm(irt/(jj) 

Walter Gautschi provided us with numerical results for the approximation of these 
two integrals with rational quadrature rules. The integrands contain infinitely many 
poles at the points p\ = u,p2 = —u,ps = 2u,p4 = —2a;, .... For ^3 these are simple 
poles but for I4 these are poles with multiplicity two. For the orthogonal quadrature 
rule we choose the parameters ti in such a way that they match up with the poles 
of the integrand, i.e., 

ti = 0, t i+ i = -l/pi, i = l,2,.... 

For the Gaussian rule we choose ti = —1/pi. The results for Gaussian rational quad- 
rature GR(n), orthogonal rational quadrature OR(n) and Gauss-Legendre quadra- 
ture GL(n) are given in Table 4 and Table 5 3 . Both rational quadrature rules 
perform very well for ^3, which is simple to understand since these rules explicitly 
take into account the poles of the integrand. The Gauss-Legendre rule does not 
take into account these singularities of the integrand and as a result it behaves 
poorly. 




3 These computations were done by W. Gautschi on the Cyber 205, with a requested error 
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UJ 11 


ClTttn} OTttn} Cn,(n\ 

\JtL\\ILi \JLY\lll VjrlJl/61 


2 5 
10 


0.46(-ll) 0.26(-8) 0.52(-5) 
0.16(-24) 0.79(-22) 0.10(-10) 


1.1 5 
10 
15 


0.21(-8) 0.26(-5) 0.21(-1) 
0.19(-22) 0.30(-16) 0.26(-3) 
0.81(-25) 0.67(-25) 0.32(-5) 


1.01 5 
10 
15 


0.43(-8) 0.93(-5) 0.24 
0.91(-22) 0.24(-15) 0.64(-l) 
0.40(-26) 0.61(-26) 0.17(-1) 



Table 4: relative errors for J 3 



When the integrand has poles of second order, then there is also a difference in 
the performance of the rational quadrature rules. 



u n 


GR(n) OR(n) GL(n) 


2 5 
10 
15 


0.65(-5) 0.79(-8) 0.47(-4) 
0.12(-10) 0.24(-21) 0.19(-9) 
0.22(-16) 0.15(-24) 0.56(-15) 


1.1 5 
10 
15 


0.17(-1) 0.27(-5) 0.13 
0.32(-3) 0.30(-16) 0.34(-2) 
0.38(-5) 0.27(-24) 0.62(-4) 


1.01 5 
10 
15 


0.37 0.16(-5) 0.78 
0.93(-l) 0.40(-16) 0.40 
0.23(-l) 0.11(-24) 0.16 



Table 5: relative errors for I4 



For the integral I4 the orthogonal quadrature rule OR(n) works very well, even 
for uj close to one. Both the Gauss-Legendre rule GL(n) and the rational Gaussian 
rule GR(n) have difficulties when to approaches 1, and even for uj = 2 they do not 
converge as fast as the orthogonal quadrature rule OR(n). For GL(n) the reason is 
clear: no poles of the integrand are being taken into account. On the other hand 
GR(n) approximates the integrand by a rational function with simple poles, whereas 
the integrand actually has double poles. Under estimating the order of the poles 
thus seems to degrade the convergence. By choosing q2n-i(x) = K n (x)ir n+ i(x) / (1 + 
t k x) 2 in (2.4) we see that the orthogonal quadrature rule integrates (1 + tkx) 2 
(k = 2, 3, . . . , n) exactly. Hence the orthogonal rule is more natural when the 
integrand has poles of order two, but in this case one could also use the Gaussian 
quadrature rule corresponding to the basis function 

1111 11 
1 + hx 1 (l + tiz) 2 ' l + t 2 x' (l + t 2 x) 2 "" ' 1 + tnx' (l + t n x) 2 ' 
which by Theorem 7 uses the zeros of the orthogonal polynomials of degree n with 
weight function nr=i(l + tix)~ 2 . 
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